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Abstract 

The Dirac-Moshinsky oscillator is an elegant example of an exactly solv- 
able quantum relativistic model that under certain circumstances can be 
mapped onto the Jaynes-Cummings model in quantum optics. In this work 
we show, how to do this in detail. Then we extend it by considering its cou- 
pling with an external (isospin) field and find the conditions that maintain 
solvability. We use this extended system to explore entanglement in rela- 
tivistic systems and then identify its quantum optical analog: two different 
atoms interacting with an electromagnetic mode. We show different aspects 
of entanglement which gain relevance in this last system, which can be used 
to emulate the former. 
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1 Introduction 

The fact that the Dirac equation is analytically solvable both for the free parti- 
cle and the Coulomb problem made it plausible to look for other solvable Dirac 
systems whose non-relativistic counterparts, in analogy to the above, are super- 
integrable and/or algebraically solvable. Indeed additional problems are solvable 
in the relativistic setting; in particular, a generalization of the harmonic oscilla- 
tor, with its large dynamical symmetry pQ. Moshinsky and Szczepaniak [2] did 
not only realize this, but constructed explicit solutions. This line was extensively 
followed by Moshinsky and his collaborators [2j |3] . 



We shall show that this solvability results from an additional conserved quan- 
tum number. This conservation indeed allows to preserve solvability in an alge- 
braic sense even as we couple an additional isospin field into the system. Ap- 
propriate coupling conserves the quantum number and thus the symmetry of the 
extended problem. This solution was clarified by a mapping onto a more general 
problem in quantum optics, which we solved recently 

The connection to the quantum optics problem not only makes the solvability 
more understandable in a wider setting but also allows to conceive quantum optics 
experiments that emulate this system. Nowadays these experiments are feasible 
with trapped ions. The free particle Dirac equation has already been emulated 
[6] and the mapping of the standard Dirac-Moshinsky oscillator (DMO) in 
2 + 1 dimensions on the Jaynes-Cummings model (JCM) has been proposed [7]. 
Here we shall review these methods and also show some cases in which the 3 + 1 
dimensional example can be mapped. We then propose a system of two atoms 
interacting with a cavity with different couplings on which the DMO coupled to 
an isospin field can in principle be realized. 

The paper is organized as follows. In the section [2] we review and describe the 
Dirac-Moshinsky oscillator, in section [3] we present the Jaynes-Cummings model 
which is of great importance in quantum optics and cavity QED. In section [4] we 
describe how to map the relativistic model of the DMO onto the JCM for the 
three dimensional cases, this will allow us to formulate a straightforward method 
to extend the model as presented in section [5j In section [6] the quantum optical 
model related to the extended model will be presented. 



2 The Dirac-Moshinsky oscillator 

The DMO, introduced in 1989 by Marcos Moshinsky and A. Szczepaniak is a solv- 
able quantum relativistic model which in the non-relativistic limit corresponds to 
the harmonic oscillator plus a spin-orbit coupling term. Noting that the mo- 
mentum operator is linear in the Dirac equation, they proposed a linear term in 
position as well. We chose to write it as follows: 

= ( CC * ' ^ P + imuj P r } + m( ?P) \^) = H\^) . (1) 

Here c denotes the speed of light, m and uj are the mass and frequency of the 
oscillator, and we make use of the following Dirac matrices: 

OL = -<j y ®s, (3 = a z ®I 2 . (2) 

We use two sets of Pauli matrices, namely s = (s x ,s y ,s z ) for the spin of the 
particle and a = (a x ,a y ,a z ) for the isospin. The representation of both a and 
s is the same, but we use different symbols to avoid confusion. We also use the 
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convention of writing the isospin terms always at the left and in the following we 
shall omit to write the outer product (g) explicitly to simplify notation. 

Using the raising and lowering operators a± = \{<J X ± icr y ) and the isospin 
states |±) = a± |=f) the state vector can be expressed as 

I*) = I-) Wt + 1+) hfe) (3) 

and the Hamiltonian can be rewritten as 

H = mc 2 a z + cd-S ■ {mujr + ip) + cu+s • (mcor — ip) . (4) 

This form shows clearly the coupling between and \ip2), the big and small 
component of the state vector respectively. Squaring the Hamiltonian results in 

2 IV'i) = {p 2 + rn 2 uj 2 r 2 — 3hujmc 2 — 2mc 2 ujs ■ L) , (5) 

where we have introduced the angular momentum operator L = r x p. An 
analogous expression can be obtained for \ip2)- At this point one can take the non 
relativistic limit noting that E = mc 2 + e, with e ^ mc 2 being the non-relativistic 
energy. The term at the left part becomes approximately 2mc 2 e, which means 
that e is an eigenvalue of the operator at the right hand side which we recognize 
easily as the Hamiltonian of an isotropic harmonic oscillator plus a constant and 
a spin-orbit coupling term. 

From equation ^ it is evident that H 2 commutes with J 2 (J = L + S), 
the total angular momentum. For H one has an additional constant of motion, 
a* ■ a + \a z that becomes clearer if one express the Hamiltonian in terms of 
the oscillator creation an annihilation operators a = \/^rjf( r + ^^j) an d = 



2h V l mujJ 

^ st _L -n I rr a . n\ 



H = mc a z + -q [a^s • a T + <j + s ■ aj , (6) 

where we have defined r/ = V2mc 2 huj. This means that the eigenfunctions of H 
can be expressed as expressed as a combination of |±), the isospin components, 
and In ± 2, 3 m j)i the eigenstates of 3D harmonic oscillator coupled to spin 
\. n is the radial quantum number, j is the angular momentum quantum number 
and rrij its projection. 

Before going into any further details of the DMO eigensystem, we take a short 
detour to introduce the quantum optical model related to our construction. 

3 The Jaynes-Cummings model 

The JCM [8] is a paradigm in quantum optics which can be thought of as a 
model that describes the interaction between a two level atom and one mode of 
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the electromagnetic field, i. e. a cavity mode. The Hamiltonian describing such 
system can be written as 



Hj C = tt(a+a + a+a)) + 5a z , (7) 

where at and a represent the creation and annihilation operator of photons in 
the cavity and a± are the raising and lowering operators of the atomic states. f2 
represents the coupling strength between the cavity and the atom and S stands 
for the detuning of the atomic transition frequency from the cavity mode. The 
energy term of the mode is absent from ([7]) , because it is written in the interaction 
picture. 

To diagonalize the Hamiltonian ([7]) one notes that there is a conserved quantity 
I = a^a+^CTz which can be identified as the number of excitations in the system as 
there is only excitation exchange in the system. Using the basis |— ) |n), |+) |n — 1), 
where |— ), |+) are the ground and excited state of the atom and |n) is a number 
state in the cavity, the Hamiltonian is block-diagonal in terms of 2 x 2 matrices: 



#ic(n) = 




The eigenenergies can be easily obtained: 

£±(n) = ±VS 2 + n 2 n = ±£(n) (9) 

and one can find the corresponding eigenstates, which are known as dressed states 
in the literature: 



\ip+(n)} = sin (0„) |-) |n> + cos (6 n ) |+) |n- 1) 

|¥>_(n)) =cos(0 n )|-}|n}-sin(0 n )|+}|n-l) (10) 



with 



Now that we have fixed the notation for the JCM and shown the simplicity of its 
solutions, we proceed to relate it with the DMO in the next section. 



4 Mapping the DMO onto the JCM 

In this section we describe the connection between the DMO and the JCM and 
show under which circumstances the DMO can be mapped onto the JCM. 
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4.1 1 + 1 DMO 



Let us consider only one spatial dimension, namely the 1 + 1 DMO. For this case 
one needs only two anticommuting Dirac matrices and we choose to write 

= —ca y (p + imu)a z x) + mc 2 a z . (12) 

where the superscript in the Hamiltonian indicates that we are working in the one- 
dimensional case. Using the creation and annihilation operators a\ = \/llF x — 
*mxJ' ax = \/ T1 W x ~^mb an d the raising and lowering operators a± = (a x ±a y )/2 
of the Dirac Spinor one can rewrite the previous equation as 

= V2mc 2 huj (a + a x + <7-a\\ + mc 2 a z . (13) 

If one takes a look at equation 0, the connection in this case is obvious. This 
Hamiltonian is exactly the JCM Hamiltonian in quantum optics. Thus the 1 + 1 
DMO maps exactly onto the JCM, provided that one identifies V2mc 2 huj — > f2, 
mc 2 — > 8, the isospin with the atomic system and the spatial degrees of freedom 
with the cavity mode. 



4.2 2 + 1 DMO 

Now let us consider the case in two spatial dimensions, the 2 + 1 DMO. For a 
full description of this case see [9] . Here we need three anticommuting Dirac 
matrices and we choose 

= —ca x (p y + imc 2 a z y) — ca y (p x + imc 2 a z x) + mc 2 a z (14) 

The ladder operators for each spatial dimension x and y can be used to construct 
a chiral representation in the form 

(a x -ia y )/\/2 a] = (a x + ia\)j\p2 



a r = (a x + ia y ) I W2 a\ = (a x — ia y )/V2 (15) 

which are also creation and annihilation operators with the canonical commuta- 
tion rule [a r ,aj] = 1 and [a/,a[] = 1. Using the chiral ladder operators, together 
with the previously defined a±, one finds 

if (2) = 2Vmc 2 huj (a + a r + a + a\) + mc 2 a z . (16) 



The set of operators a/, a\ is absent from equation (16), which means that the 
eigenstates of depend only on the number states \n r ), of the number operator 
aja r , and are infinitely degenerate in the subspace spanned by The connec- 
tion to the JCM is also obvious. One has to identify the subspace of a r with the 
cavity mode, the isospin with the atomic system and in this case mc 2 — > 6 and 
2\/ mc 2 huj —> 
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4.3 3+1 DMO 

Returning to the 3 + 1 case described by the Hamiltonian in equation ([6]), one 
recognizes that 1^ = a) ■ a + ^cr z is a conserved quantity. This tells us that a 
natural way of labeling the eigenstates of the total angular momentum and the 3D 
harmonic oscillator, is in terms of the oscillator quantum number N = 2n + j rb ^ ■ 
We define 

\N) = \n(j±\,\)jm 3 ). (17) 

There is no ambiguity in the previous definition if one remembers the dependence 
of N on both j and n, the total angular momentum and radial oscillator quantum 
numbers. The key here is the parity of N as there are two orthogonal states with 



the same value of j and n , in eq. (17), each one of these will be labeled by an N 
with the same parity as j ± |. The lowest possible value will be iV m i n = j — |. 

In order to find the eigenstates of H one has to know how the ladder operator 
s ■ a, and its hermitian conjugates, act on the number states we just defined. We 
take the result from reference [I] and write it as 

s-a\N) = y/n(N) \N - 1) 
s • at \N - 1) = ^fx(N) \N) , (18) 

with 

f J2n + 2j + 2 N = 2n + j + \ 
KN) = { ' 2 (19) 

[ V2n N = 2n + j-\. 

which takes into account the two separate cases, when N has the parity of j + \ 
and j — \- Using the basis where 1^ is diagonal, namely |— ) \N) and |+) |iV — 1), 
one can diagonalize the Hamiltonian in equation ^ in terms of 2 x 2 matrices 

mc 2 ri\/ u(N) , 
H(N)=\ /V ^ ; I- <2U> 

r/y /J,{N) —mc 



In complete analogy with the JCM (see equation (|8j)) one can express the corre- 
sponding eigenergies in the following form: 

E±(N) = ± v / m 2 c 4 + r] 2 fj,(N) = ±E(N). (21) 

The eigenstates will have the same functional form as the dressed states of the 



JCM in equation (10), but with the number states \N) and energies E(N). 

The eigenstates with N = 2n — j + ^ are infinitely degenerate as they have 



energies that do not depend on j as can be seen from (21). In addition, the 



blocks in equation ( 20 ) have the same form to those of the JCM in equation 
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Therefore we can state that the infinitely degenerate part of the 3 + 1 DMO can 
be mapped to the JCM, if one identifies y2r] —¥ f2 and mc 2 — > 5. 

If one takes N = 2n+j + \ the degeneracy is finite. In this case the functional 



dependence of the blocks in equation (20) with ./V is different to the JCM, so it 



can not be fully mapped. Even so, one could still emulate the finite degenerate 
part of the DMO in a JCM if one restricts only to one 2 x 2 block of H( N \ here 
one should identify r/ v / 2n + 2j + 2 with fi-^/n. 

5 A DMO coupled to an external field 

Now we present the extension to the DMO interacting with an isospin field mod- 
eled as a potential that is summed to the total Hamiltonian 

H = H + §. (22) 

Among the many choices which preserve the integrability of the system, we use 
the simplest (i.e. linear) one by way of example, namely 

<$> = x (a'_A +a' + A) + W z , (23) 

Here A represents the ladder operator for each dimensionality we considered in 
the previous section, and we have denoted with primes the isospin operator that 
acts on the field degrees of freedom. The full Hamiltonian is given by 

H = rj{a-A ] + a + A) + x(°'-A ] + a' + A) + mc 2 a z + ja' z . (24) 

For a physical discussion on and the covariant form of this system see |10j . 
We shall proceed from here taking A as any ladder operator, that satisfies 

A\Af) = W)\M-l). (25) 

Next we note that, due to the additional isospin, one has the integral of motion^ 

/ = jtA + \ (a z + a' z ) . (26) 

Using the basis where I is diagonal, namely 

|-)|-')|AT+1) \+)\-')\M) |-> |+') |A0 |+)|+'> IAA-1), (27) 

the Hamiltonian is now block diagonal with its blocks given by the 4x4 matrices 

/ -mc 2 - 7 x/CAf+l) rjf(M+l) \ 

X/CAT+l) l-mc 2 V f(Af) 

77/CAf + l) mc 2 - 7 x/CAO 

\ rjf(Af) Xf(ff) mc 2 + j J 



H{N) 



(28) 



x It will be further defined for each dimensionality considered 
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where f(J\f) depends on the dimensionality one choses. 

The analysis in the previous section has allowed us to construct the general- 
ization of DMO oscillator, in any of the three dimensionalities considered, coupled 
to an external isospin field. Table [5] shows for each dimensionality the correspon- 
dence of the ladder operator A, the integral of motion /, the quantum number J\f 
and the function / introduced in equation (28). 



Table 1: Correspondence for each dimensionality. 



DMO 


A 


Conserved quantum number / 


/(AO 


W) 


1 + 1 


a x 


ala x + \{a z + o' z ) 




\n) 


2 + 1 




ala r + \(a z + a' z ) 




K) 


3 + 1 


s ■ a 


at ■ a + \{a z + a' z ) 


\M*0 


\N) 



The systems are again integrable and one can find the eigenenergies by diag- 
onalizing each block H{N). We shall not write the general solutions here as they 
can be found in [10]. Instead consider the evolution of a simple initial state and 
evaluate the entanglement of the DMO with the field. 



5.1 Entanglement with the field 

In this section we analyze the dynamical features of a Dirac particle under the 
influence of the external field. To this end and for simplicity, we use a product 
initial state formed by the lowest eigenstate of the DMO times the upper state of 
the field, namely 

l^o) = |-)|+ / >|0) (29) 



With our choice N = it follows from equation (27) that the basis is reduced 
to three states, because M has to be non- negative. Thus the evolution will stay 
confined in a 3 dimensional subspacd^ of the entire Hilbert space and the state 
vector at any time can be written as 

= B 1 (t) I-) |-'> |1) + B 2 (t) |+) |-'> |0) + B 3 (t) I-) |+'> |0> . (30) 

If one simplifies things even more, by setting r\ = % and mc 2 = 7 one can find 
simple explicit solutions for these coefficients, namely 

Bi{t) =/ (t) 

Bi(t) =\ (l - hit) + (-l) l g(t)) 1 = 2,3 (31) 



2 For J\f > one would have 4 dimensional subspaces. 
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Figure 1: Purity of the field as a function of time and the rest energy 7 = mc 2 , 
with an initial state formed as a product of the lowest DMO eigenstate E = 
and the upper state of the isospin field |+). X = V- The period of oscillations 
increase with 7 an effect already found in [10] which holds only for the resonance 
7 = mc 2 . 



with the definitions 
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fo(t) =-^sm(tj) 

T 

9{t) = 1 ± 1 cos (t (7 - 7)) + 1 cos (t (7 + 7))- (32) 
2 7 7(7 + 7) 

Our next task is to find the entanglement with the field. To this end one has to 
evaluate the reduced density matrix of the field by taking a partial trace over the 
DMO degrees of freedom. One finds the density matrix 

, |-B 1 mi 2 + \B 3 (t)\ 2 , 
p'(t) = Tr DM o{|^)) <** }= _ (33) 

\B 2 (t)\ 2 
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The entanglement of the DMO with the field can be measured by the purity [11] . 
obtained as 

P ¥ = Tr{p' 2 } = \ + \(g{t)-f {t)) 2 (34) 

Figure [l] shows the purity of the field as a function of t and 7 = mc 2 , the field 
strength equal to the rest mass energy in this case. The purity starts at value 
Pf = 1 as expected for an initial product state. Full entanglement of the isospin 
with the D occurs when purity reaches it minimum value of Pf = 1/2. Note that 
the oscillations between fully entangled and totally pure situations have a period 
which increases with 7. Normally one could expect an increase of the frequency 
with 7, but as already observed in [10], this is an effect of the resonance 7 = mc 2 . 



6 Two atoms inside a cavity 

In this section we explain the connection of the coupled DMO of section [5] with 



quantum optical systems. The Hamiltonian in equation (24) can also be used in 
the context of quantum optics to describe a system composed of two atoms inside 
a cavity, provided one identifies A with a the ladder operator of the cavity and 
each isospin with an atom. Besides one has to regard r/ and \ as the coupling of 
each atom to the cavity and mc 2 and 7 as the detuning of each transition level 
with the cavity mode. 

This means that the extended model can also be mapped in cavity QED. Again 
the 1 + 1 and 2 + 1 cases can be mapped exactly, whereas the 3 + 1 case can only 



be reproduced for the special situation considered in subsection |5.1[ because in 
3 + 1 dimensions the coupling with the field mixes the dynamics of the infinitely 
degenerate part with the finitely degenerate part. This case f(N) in the blocks 



(28) of H to depend on fi(N) which changes by steps of 2j + 2 as N increases. 



For N = the equation (28) reduces to a 3 matrix and only the value /x(0) = 2n 



comes into play. Thus the particular case studied in subsection 5.1, which can 
represent the three dimensionalities, can be mapped to a quantum optical system. 
For this reason, in this work we focus only on the equivalent quantum optics model 
restricted to n = 0, that is an initially empty cavity. 

In the context of quantum information theory, this is one of the simplest mod- 
els which can be used in order to study two important aspects: entanglement 
and decoherence. The former is a resource for implementing quantum informa- 
tion protocols - and therefore, our ally - while the latter is an obstacle for such 
implementation. While these concepts play an important role in nurturing the 
quantumness of a central system, here we shall make use of quantum informa- 
tion measures as a tool to infer dynamical features of our problem in the simplest 
possible way. 



We shall not rewrite the Hamiltonian as it is equal to the one in equation ( 24 ) 



and the same solutions apply. The initial state to consider here has to be different 
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as the quantities of interest are others. We distinguish between the two atoms 
(central system) and the cavity (environment). So we shall begin with a product 
state of an atomic state times a number state of the cavity, and for simplicity we 
again take M = 



(cos (a) |-+) +sm(a) |+-)) |0) . 



(35) 



note that for a = the state is exactly the same as in [29j We shall consider 
equivalent conditions as in subsection |5.1[ which in this case means equal couplings 
to the cavity and atomic transition frequencies equally detuned from the cavity 
mode. With these considerations at time t the state vector can be found in in a 



superposition like in equation (30), the only difference is that the coefficients will 



depend on a, so we shall write them as 
Ci{t) =f a (t) 

Q(t) =\ (l - f a (t) + (-l) l g(t) cos (2a)) , 1 = 2,3 
where g(t) and fo(t) are given in equation ( J32| ) and we introduced 

/ Q (t) = (l + sin(2a))/ (t). 
One has to remember that here 7 represents the detuning of both atoms. 

6.1 Entanglement measures 

To evaluate the measures of entanglement one needs of the reduced density matrix 
of the two atoms, so we trace over the oscillator degrees of freedom to get: 



(36) 



(37) 



p = Tr 0sc {!*(*)> <*(i)|} 



( |Ci(t)| 2 \ 

|c 2 (t)| 2 (c(t) 3 yc 2 (t) 

(c(t) 2 yc 3 (t) \c 3 (t)\ 2 

V j 



(38) 

To measure the entanglement between the two atoms (central system) and the 
cavity (environment) we use the purity P = Tr {p 2 } and find 

P(t) = l-2f a {t) + 2fl{t) (39) 

which serves as measure of the decoherence of the two atoms system. 

To measure the entanglement between the atoms we use the concurrence |12j 

1 /2 

C(p) = Max {0, Ai — A2 — A3 — A4}, where Xj are the eigenvalues of [pa y a' y p* cr y a' y ) 
in non- increasing order. In this case we find 



C(t) = V(/a(t)-l)-S 2 (*)cos2 (2a) 



(40) 
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Figure 2: Concurrence versus purity plane for an initially empty cavity and two 
equal atoms. The black curves show the case with zero detuning 7 = and initial 
state a = tt/4 (see eq. ([35])), (a = tt/40 bottom curve on the right.) In red 
the curves for detuning 7 = 1 for a = on the left and a = 7r/40 on the right. 
The gray area corresponds to unphysical density matrices and its lower frontier 
represents the so called maximally entangled mixed states. 
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A useful way of visualizing both dynamics together is the so called CP plane 
131 -4j. Figure 6.1 shows this plane for two separate cases, both with detuning 
1. The left part shows in red for an initial state with a = 0, the equivalent 
case to subsection 5.1, whereas the right part shows it in red for a = tt/40. On 
both sides the red curves are parametrized up to t = 30. The black curve on 
both figures shows the behavior for zero detuning for a = tt/4 and also on the 
right part for a = 7r/40 (lower black curve). Both figures show a gray area that 
correspond to states with density matrices which are not physically acceptable 
and its lower frontier corresponds to the maximally entangled mixed states. 

It is interesting to note that in both cases the red curves form Lissajous like 
behavior in the CP plane. For zero detuning analytic expressions for the concur- 
rence as a function of the purity can be found 



C±(P;a) 



1 ± • N /1 + 2(P-1) - 2 sin (2a) 



(41) 



these curves form the frontier of the region that the red curve partially fill as 



one can see in the figure 6.1 The region filled by the red curves reduces as one 
increases the detuning 7, that is the atoms feel less the effects of the cavity. 

We also found similar behavior in [1] when instead we had zero detuning but 
interaction between the atoms. 
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7 Conclusions and final remarks 



The Dirac Moshinsky oscillator in 1 + 1 and 2 + 1 dimensions was mapped to 
the Jaynes-Cummings model. For the 3 + 1 case at least the infinitely degenerate 
part can be mapped. This scheme was based on explicitly using an invariant of 
these systems, which in the optical image acquires the simple meaning of the total 
number of excitations, which differs from the total energy because of detuning. 

Based on this picture we obtained a soluble extension to the DMO coupled 
to an isospin field, if this coupling is carefully chosen to have an invariant that 
includes the excitations of the isospin field. The system not only retains solvability, 
but can be mapped to 2 atoms inside a cavity again for the 1 + 1 and 2 + 1 
dimensional cases. 

In the 2+1 dimensional DMO the presence of an external field which allows 
integrability does not break the infinite degeneracy. The absence of the operator 
ni for the JC model can be thought as "inert" field modes which do not interact 
with the atoms in our system. In the quantum-optical system, other cavity modes 
may be populated by photons. However, if their frequencies are far off-resonance 
with our two- level atoms, the Hamiltonian in the interaction picture will contain 
no terms related to such modes. To complete the analogy, a Dirac particle in 
2+1 dimensions contains such degrees of freedom, but they have no effect on the 
energies of the DMO. 

We recognized that in our setting for 3 + 1 dimensions, only the block N = 
can be represented by this diatomic model. However, the additional structure in 
the relativistic model may lead to other soluble models in atom optics, known or 
unknown at present. This will have to be the subject of further investigations. 

By way of example we evaluated the purity of the field by tracing over the 
DMO degrees of freedom as well as other examples of the entanglement measures 
of the quantum optical system to which we mapped many of the discussed models. 
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